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Abstract 

This paper deals with the long-time behavior of solutions of nonlinear reaction- 
diffusion equations describing formation of morphogen gradients, the concentration 
fields of molecules acting as spatial regulators of cell differentiation in developing tis- 
sues. For the considered class of models, we establish existence of a new type of 
ultra-singular self-similar solutions. These solutions arise as limits of the solutions 
of the initial value problem with zero initial data and infinitely strong source at the 
boundary. We prove existence and uniqueness of such solutions in the suitable weighted 
energy spaces. Moreover, we prove that the obtained self-similar solutions are the long- 
time limits of the solutions of the initial value problem with zero initial data and a 
time-independent boundary source. 



1 Introduction 

In the studies of reaction-diffusion equations, one canonical problem deals with the follow- 
ing equation [2,11: 



u t = Au-u p , (x,t) G R d x (0,oo). (1) 

Here p > 1 is a constant and u = u(x, t) > can be viewed as the concentration of a 
chemical species diffusing in the ci-dimensional space subject to degradation whose rate is 
an increasing function of the species concentration. Usually, one considers the associated 
Cauchy problem with some non-negative initial data u(x, 0) = uq(x). During the 1980's, 
this problem attracted a considerable attention, in particular in the case of measure- valued 
initial data (e.g., when uq is a Dirac mass) (2}§§[13KT7K24]. In the course of these studies, 
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it was discovered that ([I]) possess self-similar solutions for all 1 < p < (2 + d)/d, which 
are smooth for all t > and converge to zero outside the origin, while blowing up at 
the origin when t — > + (3j[TT] (see also [8] for a variational approach). These solutions 
play important roles in determining the long-time behavior of the solutions of the Cauchy 
problem for general classes of initial data and in some sense describe the transient dynamics 
in systems described by ([!]) |4j[l0j[TT|[T5}[T7|[24||3l] . In particular, a special class of self- 
similar solutions of called very singular solutions attract the physically important class 
of initial data with sufficiently fast asymptotic decay [9 11, 17 



Equation ([!]) with p > 1 on domains with boundaries also arises as a canonical model of 
morphogen gradient formation (for recent reviews, see |20p5p7{|30] ). Morphogen gradients 
are concentration fields of molecules acting as spatial regulators of cell differentiation in 
developing tissues [22]. In particular, the case p > 1 was proposed to describe a robust 
patterning mechanism whereby morphogen increases the production of molecules which, 
in turn, increase the rate of morphogen degradation |7|. For example, a protein called 
Sonic hedgehog (Shh) is known to induce the expression of its receptor Patched, which 
both transduces the Shh signal and mediates Shh degradation by cells in the Drosophila 
embryo [5][l6]. 

An important aspect of morphogen dynamics is the presence of localized sources at the 
boundary of the morphogenetic field. This leads to the need to consider initial boundary 
value problems, whose prototype is the following one-dimensional problem: 

u t = u xx -u p (x,t) £ [0, oo) x (0, oo), 

u x (0,t) = -a tG(0,oo), (2) 
u(x,0) =0 x £ [0,oo). 

This problem can be viewed as an extension of the Cauchy problem for defined for x > 
in the presence of a boundary source at x = 0. Here a > is a constant characterizing the 
source strength of morphogen production, and the zero initial condition corresponds to the 
absence of the morphogen at the onset of patterning. In what follows, we will restrict our 
attention only to this simplest model of morphogen gradient formation. 

In the context of morphogenesis, one is often interested in the establishment of a sta- 
tionary morphogen profile and the transient dynamics that leads to it. The stationary 
problem for Q can be written as the following boundary value problem: 

v X x -v p = 0, v x (Q) = -a, v(oo) = 0, (3) 

whose unique solution for any p > 1 is explicitly given by 

, , /2(p + l)\?=i. , ?_ /2i+T(p + l)i+l\ _P-L 

< x ) = { (^3ipJ (a + x) p-i, a=\ ^—J—\ a ^i. (4) 

In fact, it is easy to see that the stationary solution v(x) in ^ is the limit of the solution 
u(x, t) of ([2]) as t — > 00 for each x > 0, and is approached monotonically from below [14] , 
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Figure 1: Numerical solution of ^ in self-similar variables for p = 2 and a = 1. Thin lines 
show snapshots of the solution corresponding to t = 0.1, 1, 10, 100 (the direction of time 
increase is indicated by the arrow). The bold line shows the asymptotic solution. 



However, as we noted in [14] , this approach is not uniform in x and for each fixed x > 
occurs on the diffusive time scale r p (x) = 0(x 2 ), which diverges as x — > oo. Thus, the 
timing of the establishment of the steady morphogen concentration at a given point depends 
rather sensitively on the location of that point. 

To better understand the dynamics of the approach of the solution of ^ to the sta- 
tionary solution, we undertook numerical studies of the initial boundary problem in ^ 
for various values of p > 1. In those studies, we discovered that when the ratio of the 
solution at a given x to the value of the stationary solution at x is plotted vs. the diffusion 
similarity variable x/\/t, the numerical solution approaches some universal limit curve de- 
pending only on the value of p (23) . This process is illustrated in Fig. [TJ where the results 
are presented for the biophysically important case p = 2. This observation suggested to 
us some hidden self-similarity in the behavior of solutions of ^ Jl]. A simple scaling 
argument indicates that the long-time behavior of the solution of (pi) for a fixed value of 
a > is closely related to the behavior of solutions of ^ at fixed x > and t > as 
a —> oo |23|. We found numerically that in the limit a — > oo the solutions of ^ attain 
a self-similar profile (see the following section for precise definitions) [23]. The purpose 
of this paper is to substantiate these numerical observations by establishing existence and 
properties of what we will call ultra-singular self-similar solutions in the limit of infinite 
boundary source strength. We also prove that these solutions are indeed the long-time 
limits of the solutions of ^ in the above sense. 

We note that the solutions constructed by us form a new class of self-similar solutions 
to Q in d = 1. Indeed, our solutions can be trivially extended to the whole real line by 
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a reflection and can be viewed as singular solutions of ([I]) that blow up at the origin. We 
point out that these solutions are different from the self-similar solutions studied in (3 11 . 
The ultra-singular solutions of ([I]) constructed by us can be viewed as the more singular 
counterparts of the very singular solutions of [3] in the following sense: the singularity in the 
former is concentrated on a half-line (x = 0, t > 0) in the (x, t) plane, while the singularity 
in the latter occurs only at a single point {x = 0, t = 0). Similarly, our convergence result 
for the solutions of ^ with a £ (0, oo) may be viewed as a counterpart of the result of [l~7] , 
in the sense that in the former case the solution can be viewed as the distributional solution 
of ([!]) with an added term 2ad(x) in the right-hand side, while in the latter case one can 
think of the solution as the distributional solution of ([I]) with the term a5(x)5(t) added to 
the right-hand side. 

Before concluding this section, let us briefly mention a few possible extensions and open 
problems related to our present work. It would be interesting to understand the role our 
self-similar solutions play for the singular solutions of the initial value problem associated 
with ([I]) for general non-zero initial data. Let us point out that even the basic questions of 
existence and uniqueness of such singular solutions for the considered parabolic problems 



in suitable function classes are currently open (see 29 for a very recent related work). 
Other natural extensions include higher dimensional versions of the considered problem, 
as well as a proof of global stability of self-similar solutions. These studies are currently 
ongoing. From the point of view of applications, it is also important to consider solutions 
of ([I]) with added time- varying singular sources, for which both the very singular and the 
ultra-singular solutions may be relevant. 

Our paper is organized as follows. In Sec. [2j we introduce a singular version of the 
initial boundary value problem in ^ and prove existence, uniqueness, monotonicity and 
limiting behavior of the self-similar solution to this singular problem. Then, in Sec. [3] we 
prove that the obtained self-similar solutions are the long-time limits of the solutions of 
([2]) in an appropriate sense. 



2 Singular solutions and the similarity ansatz 

Let us consider ([2]) with infinite source at the boundary, i.e., the following singular initial 
boundary value problem: 

u t = u xx -u p (x,t) G (0,oo) x (0,oo), 

u(0,t)=oo, te(0,oo), (5) 
u(x,0) = x € (0,oo). 

By a solution to ([5]), we mean a classical solution for all (x,t) 6 (0,oo) x (0, 00) decaying 
sufficiently fast as x — > +00 for all i > 0, and continuous up to t = for all x > 0. Note 
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that for each p > 1 this problem possesses a singular stationary solution 



(p — l) 2 / \x 

which is the limit of v a (x) as a — > oo for each cc > 0. 

Consistently with the discussion in the introduction, we now seek solutions of §5§ in 
the form 

u(x, t) = uoofcMC), C = Mx/y/i), (7) 

for some function < <p(C) < 1, which will be referred to as the self-similar profile. 
Substituting the similarity ansatz from ([7]) into ([5]), after some algebra we obtain the 
following equation for the self-similar profile 0: 



which must hold for all £ G (— oo, oo), supplemented with the limit behavior 



(8) 



lim 0(C) = 1, lim ^ = 0, (9) 

lim 0(C)=O, lim ^P=0. (10) 

C->-+oo C^+oo dC 

Existence and multiplicity of solutions of ^ satisfying Q and ( 10 ) are not at all a priori 
obvious in view of both the non-linearity and the presence of singular terms in the consid- 
ered boundary value problem. In 23 , we were able to construct such solutions numerically 
for several values of p. Here we establish their existence and uniqueness for all p > 1 within 
a natural class of functions. 

We will prove existence and uniqueness of solutions of ^ satisfying ^ and (10) in 
the weighted Sobolev space ff 1 (IR, dp,), which is obtained as the completion of the family 
of smooth functions with compact support with respect to the Sobolev norm H-llffi^,^), 
defined as 

IMlW,^) = IKHi 2 (iMM) + IMIi 2 (w ( n ) 

where | \w\ |^2( K df A = f^ w2 {C)dp{(), and the measure dp is 

dp(C) = p(()d(, P (C) = exp | ^ - f |±1) c} • (12) 
Our existence and uniqueness result is given by the following theorem. 
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Theorem 1. There exists a unique weak solution (ft of Q, such that 



- r\ £ fl^d/i), 
), suc/i i/iai 77(C) = 1 /or all ( < and r/(() = 
, satisfies ([8| classically and < (ft < 1. In 



with p defined in (12), for every r] € C c 
for all Q > 1. Furthermore, (ft £ C c 
addition, (ft is strictly decreasing and satisfies ^ and (10). 

Before proceeding to the proof of Theorem [TJ let us establish a basic technical lemma 



needed to deal with the weighted spaces introduced above, which is an extension of 21 
Lemma 4.1] for exponentially weighted Sobolev spaces (cf. also (8| Lemma 1.5]). 

Lemma 1. Let w G -ff 1 (M, dp). Then there exists Rq > such that 



-2R 



w dp, < 



R 



(dwY 



dp 



Vi? > #0, 



(13) 



and 



p(R)w 2 (R) < 2e 



-R 



R 



f ) ^ 



Moreover, there exists R' < such that 



for a.e. R> Rq. 



(14) 



R 



w 2 dp < 8 



p-1 
p + 3 



dw\' 

dc) 



dp 



Vi? < R'r 



(15) 



and 



p{R)w 2 (R) < 8 



p — 1 
p + 3 



I)* 



/or a.e. R < R' Q . 



(16) 



Proof. Arguing by approximation, observe that by an explicit computation and an 
application of Cauchy-Schwarz inequality we have 



w 2 (R)p(R) + 



R 



00 dw , 
w—rrdp < 



R 



d( 




In p J w dp 



w 2 dp 



R 



dw\ 



2 

I dp 



1/2 



(17) 



In particular, (12) and (17) yield 



,2R 



p + 3 
p-1 



w 2 dp < 4 



R 



R 



dw 
d( 



dp, 



(18) 



which for large enough R implies ( 13 ). Next, since ^ In p > for large positive £, dropping 
the second term in the left-hand side of (17) and using (13), we obtain (14). 
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Similarly, we note that 



w 2 {R)p(R) 



R 



dC 



Inp ] w 2 dp 



R dw , 
W—rzdp < 



which implies 



P + 3 
p — 1 



d( 



2R\ 2 r R 



w 2 dp 




(19) 



w 2 dp < 4 



fdw\' 



(20) 



and thus ( 15 ) holds for sufficiently large negative R. Finally since J? In p < for large 



negative £, from (1 191) and (15) we obtain (16). 



□ 



Proof of Theorem [T| The proof consists of five steps. 

Step 1. We first note that ([8| is the Euler-Lagrange equation for the energy functional 




<ft 2 (p + l - 2<pP~ l ) 



(21) 



where r/(£) is as in the statement of the theorem. Indeed, the functional £ in (21) is 
continuously differentiable in i7 1 (M, dp) in the natural admissible class A defined as: 

A := {(ft G #ioc(K) : (j) — 7] £ fl" 1 (R, <ip), < < 1}. (22) 



Note that the role of 77 in the definition of £ is to ensure that the integral in (21 ) converges 
for all (ft £ A. The precise form of r](() is unimportant. Then it is easy to see that the 
weak form of ^ in i? 1 (M, dp) is precisely the condition that the Frechet derivative of £[<$>] 
is zero. 

Step 2. We now establish weak sequential lower-semicontinuity and coercivity of the func- 
tional £ in the admissible class A in the following sense: let 4>k = V + w k, where Wk — w 
in Ji' 1 (R, dp). Then 1) liminffc^oo £[(j>k\ > £[</>], where (ft = rj + w, and 2) if £[(/>&] < M for 
some MGl, then 1 1 w?fc 1 1 ifi (r,^) < M' for some M' > 0. 

Let us introduce the notation £[</>, (a, b)] for the integral in (21), in which integration 
is over all £ G (a, 6). Then, using ( |13| ) from Lemma [T] we find that for 12 > 1 



(fl, +oo)] > e 



2B 



p + 1 

(p-l)3 



w k dp > 0. 



(23) 



Similarly, taking into account that the integrand in (21) is non- negative for £ < 0, we 
have £[(ftk, (— oo, —12)] > for every R > 0. Since £[•, (—12, 12)] is lower-semicontinuous by 



7 



standard theory [6j, we obtain £[4>k] > £[4>k, yielding the first claim by passing 

to the limit R — > oo. 



To prove coercivity, we first note that by (13) 

£[<t> h ,(R,+oo)] > I 
Jb 





dw k 




dC 


1 


/•oo 


> - 




- 4 


JR 



p+1 2 I , 



(24) 



for large positive R. On the other hand, since p — 1 — 4> 2 {p + l — 2(j> p 1 ) > (p — 1)(1 — 4>) 2 
for all < d> < 1, we have 



Finally, by boundedness of 0^ and 77, we also have 



£[ct> k ,(0,R)} > 




+ w\ > d/j, — CR, 



(25) 



(26) 



for some C > independent of Wk- So the second claim follows. 

Step 3. In view of the lower-semicontinuity and coercivity of £ proved in Step 2, by the 
direct method of calculus of variations there exists a minimizer (ft £ A of £. Noting that 
since the barriers (f> = and = 1 solve ^ as well, we also have (see e.g. |19|) that <j) 
is a weak solution of Q by continuous differentiability of £ in H 1 ^, dfi) noted in Step 
1. Furthermore, by standard elliptic regularity theory [l2j, cj) £ C°°(M) and is, in fact, a 
classical solution of Q. Also, by strong maximum principle [12] , we have < (f> < 1. To 
show monotonicity, suppose, to the contrary, that 4>(a) < cj)(b) for some a < b. Then </>(£) 
attains a local minimum for some Co £ ( — oo, 6). However, by ([8]) we have d 2 4>(Co) I d( 2 < 0, 
giving a contradiction. By the same argument difi/dQ = is also impossible for any ( S 1. 
Finally, since (j) — rj € i/ 1 (M, cZ/i), monotonicity implies the first condition in ^ and (10). 

Step 4. We now discuss the asymptotic behavior of minimizers obtained in Step 3 as 
C — > ±oo and, in particular, prove the second parts of ^ and (10) and the fact that every 
solution of Q belonging to A has the same asymptotic decay, which will be needed later. 
Let us first consider the case of £ — > +oo. Performing the Liouville transformation by 
introducing 

<j)^p-£L 2 (R,+oo), (27) 



if, 



where p is defined in ( 12 ) and R > 1 is arbitrary, we rewrite (JSj) in the form 



dC 2 



£>R. 



(28) 



S 



Here q(() = g (C) + ?i(C), wnere 



e 4 ^ » — 5 
4 p - 1 



: 2C + 1 , 



*C0 = ^ ^CO- 



CM) 

(30) 



Observe that > go(C) > 1 > f° r an ( > R, with i? sufficiently large positive. 
Therefore, (28) has two linearly-independent positive solutions ipi and ifa, such that ipx — >• 
and ^2 — ^ oo together with their derivatives as C - > +°° (see e.g. |28|). In particular, 
ip = Cipi G L 2 (R, +oo) for some < C < oo, and by direct computation 



dip C ( dipi ipt d 

dC = ^/p\~dC~YdC np 



as C ~~ ► +oo. 



On the other hand, as follows from (14), we have 

Qi (0 = o{p 1 -^), 



(31) 



(32) 



so <7i(C) nas a super-exponential decay as Q — > +oo. Let ipo be the unique positive solution 
of (28) with q = qo and ipo(R) = 1 which goes to zero as ( — > +oo. Then we claim that 
■0i(C)/V'o(C) ~~ ► c f° r some < c < oo. Indeed, functions V>i and satisfy 



dC 2 



fa>(C) + 9i(C)Wi, 



dC 2 



C > 



(33) 



Multiplying the first and the second equation of (33) by and ipi, respectively, and taking 
the difference, we obtain 



d_ 



, dipt dipo , 



(34) 



Integrating this equation and taking into account that ipo,ipi and their derivatives vanish 
as £ — 7- +oo, we have 



and therefore 



hi 



Integrating this equation again, we obtain 



In 



qi(s)ip (s)ipi{s)ds, 

c 

.ipi(s)ip (s) 
gi(*) , , , . da. 



<7i (s) . , ^ , , ^ dsda. 



R Ja 



(35) 



(36) 



(37) 
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In a view of boundedness of functions ipo and we have \jpo(s)/tpo(a)\, \ipi(s) /ipi(a)\ < C 
for some C > and all s > a > R. Moreover, the estimate in (|32[) gives |<7i(s)| < C exp(— s) 



for some C > and all s £ [R,oo). Therefore, the integral in the right-hand side of (37) 
converges: 



R J a 



91 (*) 



^i(s)^o(s) 



^l(o-)^o(cr) 



R J a 



e- s dsda < Ce~ R 



< oo, 



(38) 



which immediately implies that the ratio of ipo and t/ji approaches a finite non-zero limit 
as C — > +oo. 

We can use a similar treatment to establish the asymptotic behavior of minimizers 
when £ — > — oo. The Liouville transformation 



(i-0)VpeL 2 (-oo,i?), 



with p defined by ( 12 ) and arbitrary R < applied to (|8|) yields 

d 2 e 



dC 2 

Here r(Q = r*o(C) + ri(£), where 



ro(C) 



n(C) 



1 

4 



r(C)^, 



3p+l 
p-1 



(39) 



(40) 



+ 



P 



i e2C + ^r 

p-l 4 



=4C 



2(p + i) (4>{i -4P- 1 ) 



l 



(p-l)2 

By direct computation, note that in the limit £ 

2 



ro(C) 



1 (3p + l 
4 1 p- 1 



-oo we have 
ri(C)-Kr. 



(41) 
(42) 

(43) 



Therefore, ro(C) > f(C) > 4 > for all £ < R with sufficiently large negative, and (40) 
has two linearly-independent positive solutions 0\ and 82 such that 6\ — > and #2 — > °o 
together with their derivatives as £ — )• —00. In particular, 9 = C9\ £ L 2 (— 00, i?) for some 
< C < 00, and 



7^ VdC "Ydc 



In p 







as C 



-00. 



(44) 



On the other hand, as follows from ( 16 ) we have 



r 1 (C)=o(p- 1/2 ), 



(45) 
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so ri(£) has an exponential decay as £ — )■ — oo. Computations practically identical to those 
presented above show that the ratio of 9q (the solution of (40) with r = tq which decays 
as C — ► — oo) and 6\ tends to a positive constant as Q — > — oo. 

Step 5. We now prove uniqueness of the obtained solution, taking advantage of a sort of 
convexity of £ similar to the one pointed out in [18] . Suppose, to the contrary, that there 
are two functions <frx,(j>2 G A which solve ((8|) . Define 



m+ii-t)^. (46) 

We claim that 0* G A as well. Indeed, in view of the result of Step 4 we have m < 0i/02 < 
M for some M > m > and, therefore, 

ll<^llL2((0,l),^) ^ C i H^llL2((l,oo),d M ) ^ ll^lllL2((l,oo),d M ) + \\M\h((l,oo),dfl)' (47) 



kt\\2 



i - m - (i 




^((-oc.O)^) J_ oo y i+ v / t0 2 + (1 _ f 

< C(\\l - 4>l\\L 2 ({-oo,0),diM) + II 1 - 02||l 2 ((-oo,O),<^)) 2 ' (48) 

< C(|l*i/rfClli"(E,« + H*/<iCllL'(E,*)) 2 , (19) 

for some C > 0. In fact, it is easy to see that the function E(t) := £[0*] is twice continuously 
differentiable for all t 6 [0,1]. A direct computation yields 



d 2 E(t) 

~ k x m\ + (i - t)<p\Y v z d( ri d( 

i? - 2 .) 2 (t0 2 . + (1 - ^i)^ Um(C). (50) 



P + 1 / -2 i2\2(.<2 



2p-2 

Therefore, d 2 E{t)/dt 2 > for all t £ [0,1], and so E(t) is strictly convex. However, since 
the map t i— > <p t — n is of class C 1 ([0, 1]; i/ 1 (IR, d/i)), which can be seen by a computation 
analogous to the one in (50), this contradicts the fact that dE(0)/dt = dE(l)/dt = by 
the assumption that 0i and 02 solve weakly Q and hence are critical points of £. □ 

Remark 1. Results of Step 4 °f the proof above allow to obtain the precise asymptotic 
behavior of the solution of (|8j) constructed in Theorem^by using the exact solutions of the 
associated linearizations o/ ([8]) about = and = 1. These asymptotics read 23h 

( e 2C 5 _ v \ 
0(C)~exp(^-— + — Cj , C^+oo, 

1 - 0(C) ~ exp (fcilc) , C^-oo. (51) 
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3 Long time behavior of solutions for problem (12 



In this section we prove that the ultra-singular solutions constructed in Sec. [2] have a 
direct relevance to the long time behavior of solutions for the problem in Q . Specifically, 
solutions of ([2]) converge to self-similar profile <j) at the fixed ratio x/\/i as t — > 00. That 
is, the following result holds: 

Theorem 2. Given a > 0, let u and v be the solutions of ^ and respectively, and 
set 



Then 



Moreover, 



lim F((,t) = </>(() VCGM. (53) 

t— >oo 



4>(0 <F((,t)< 0(C) (54) 
where £,(C,t) = hi(e^ + bt^ 1 ^ 2 ) and b is some large enough constant. 

Proof. The proof relies on a direct application of the comparison principle. We start 
with a formulation of the comparison principle which will be applied to ([2]). Define the 
following quantities 

P[u] =u t - u xx + u p , (55) 

Q[u] = u x + a, (56) 

assume that the functions u and u satisfy the differential inequalities 

P[u] > 0, t > 0, x > 0, (57) 

Q[u] < 0, t > 0, x = 0, (58) 

and 

P[u] < 0, t > 0, z > 0, (59) 

Q[u] > 0, i > 0, x = 0. (60) 

and, in addition, assume that n(x, t = 0) = n(x, t = 0) = 0. Such functions are called 
super- and sub-solutions for ([2]) and have the property [26] : 

m(x, t) < u(x, t) < u(x, t), (x,t) G [0,oo) x [0,oo). (61) 
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In what follows we will explicitly construct sub- and super-solutions for ([2]). 
We first show that the function 



u(x,t) = v(x)4>(z), z = In 



x + b 



(62) 



is a sub-solution, provided b > a is large enough. Here 4> verifies ([8]), Q and (10), and a 
is defined in Q. Direct substitution of (62) into (59) gives: 



P[u] 



v{x) 
(x + b) 2 

\ P -i 



1 



x + b 



d 



+ 



2(p + l) 



1 



x + aj \ dz^ J ' (p — l) 2 
In view of the fact that d<j)/dz < we have 

F[u] < \/x > 0, Vt > 0, 



x + 6 



x + a 



[l-tT 1 ) • (63) 



(64) 



provided that b > a. 

Next, direct computations also give 



Q[u(x = 0,t)] 



2,4 



p+i 

(p — l)ap- 1 



. . 1 a(p — 1) d , , 



(65) 



Let us show that Q[u(x = 0,t)] > for t > when 6 is large. To do so, it is enough to 
show that 



g(z) := 1 - <j)(z) + £~4>{z) > V2GI, 



(66) 



for e > small. Indeed, observe first that lim 2 ^+oo 5(2;) = 1 and lim 2 _ ! ._ OC) g(z) = 0. So, if 
(66) is violated, g(z) has a local minimum at some point z* £ R with 5(2;*) < 0. Since z* 
is a critical point we have 







d 

dz 



9{ z 1 = -l^) + e^ { z*) = 
1 + e 6 ^- - e P -^) ~4>(z*) - 2e {p + l \ (^{z*) - <f?(z*)). (67) 



2 P-IJ dz TV ' ' (p- l) 2 
Therefore, there exists e £ (0, 1) such that 

|0,(**)|<WCO, 

Thus, from the definition of g we have 

g(z*)>(l-e)(l-<P(z*))>0, 



(68) 
(69) 
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contradicting our assumption about g(z*). Finally, choosing b = max{a, 



a(p-l) 
2e 



} we have 



that the conditions in ( 59 ) and ( 60 ) are satisfied and thus ( 62 ) is indeed a sub-solution for 
it. 

Now we turn to the construction of a super-solution, which we will seek in the form 



u(x,t) = v(x)4>((), ( = ln 



Vt 



Straightforward computations give 

v(x) r 4 



P[u(x,t)} 



x 2 yp — 1 
2(p + l) 



x 



+ 



x + a 

x 



d 



(p-iy 



x + a 



0(1 -^ P - 1 )}: 



and 



Q[u(x = 0,t)} 



(p 



¥+T, hm 



, 1 a(p — 1) _ r d , , 



It is clear that P[u(x, t)] > for all t > and x > 0. Let us now show that 

Q[u(x = 0,t)] = Vi>0. 

Since by (|9j) and ^ 



lim (1 - 0(C)) = 3 

(,— > — oo 



we only need to show that 



lim e C -^<XC) 



0. 



Indeed, multiplying pi) by p we have 



d 



d 



d( \ p dC ( 



2(p+i; 



p^i-^- 1 ). 



(p-l)2' 

Integrating this equation and rearranging terms involving p, we obtain 



exp 



e 2 f 

— j 1 

4 p 



1 



^(^(l-^- 1 ^))^ 



(70) 



(71) 

(72) 

(73) 
(74) 

(75) 
(76) 



(77) 



By (51) we have p(C)<ft(()(l — (ff ~ exp(£) as £ — >• — oo and thus the integral in the 

right-hand side of (|77|) converges as £ — )• — oo, which readily implies (75). Therefore, both 



conditions (57) and (58) are satisfied and so (70) is a super-solution. 



Finally, the statement of the theorem follows from (61), (62) and (70). 



□ 
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Remark 2. Note that the result of Theorem^ may be extended to problem ([2]) in which 
the constant a is replaced by a bounded, monotonically increasing function a(t) > 0. 

Acknowledgements. This work was supported, in part, by NSF via grant DMS-1119724. 
CBM would also like to acknowledge partial support by NSF via grant DMS-0908279. We 
wish to thank S. Shvartsman for suggesting this problem to us and V. Moroz for helpful 
comments. PVG also would like to acknowledge valuable discussions with S. Kamin. 



References 

[1] G. I. Barenblatt. Scaling, self- similarity, and intermediate asymptotics. Cambridge 
University Press, 1996. 

[2] H. Brezis and A. Friedman. Nonlinear parabolic equations involving measures as initial 
conditions. J. Math. Pures Appl, 62:73-97, 1983. 

[3] H. Brezis, L. A. Peletier, and D. Terman. A very singular solution of the heat equation 
with absorption. Arch. Rational Mech. Anal., 95:185-209, 1986. 

[4] J. Bricmont and A. Kupiainen. Stable non-Gaussian diffusive profiles. Nonlinear 
Anal, 26:583-593, 1996. 

[5] Y. Chen and G. Struhl. Dual roles for patched in sequestering and transducing hedge- 
hog. Cell, 87:553-563, 1996. 

[6] G. Dal Maso. An Introduction to T -Convergence. Birkhauser, Boston, 1993. 

[7] A. Eldar, D. Rosin, B. Z. Shilo, and N. Barkai. Self-enhanced ligand degradation 
underlies robustness of morphogen gradients. Devel. Cell, 5:635-646, 2003. 

[8] M. Escobedo and O. Kavian. Variational problems related to self-similar solutions of 
the heat equation. Nonlinear Anal., 11:1103-1133, 1987. 

[9] M. Escobedo and O. Kavian. Asymptotic behaviour of positive solutions of a nonlinear 
heat equation. Houston J. Math., 14:39-50, 1988. 

[10] M. Escobedo, O. Kavian, and H. Matano. Large time behavior of solutions of a 
dissipative semilinear heat equation. Comm. Partial Differential Equations, 20:1427- 
1452, 1995. 

[11] V. A. Galaktionov, S. P. Kurdyumov, and A. A. Samarskii. Asymptotic "eigenfunc- 
tions" of the Cauchy problem for a nonlinear parabolic equation. Mat. Sb. (N.S.), 
126:435-472, 1985. 



15 



[12] D. Gilbarg and N. S. Tmdinger. Elliptic Partial Differential Equations of Second 
Order. Springer- Verlag, Berlin, 1983. 

[13] A. Gmira and L. Veron. Large time behaviour of the solutions of a semilinear parabolic 
equation in ~R N . J. Differential Equations, 53:258-276, 1984. 

[14] P. V. Gordon, C. Sample, A. M. Berezhkovskii, C. B. Muratov, and S. Y. Shvartsman. 
Local kinetics of morphogen gradients. Proc. Natl. Acad. Sci. US., 108:6157-6162, 
2011. 

[15] L. Herraiz. Asymptotic behaviour of solutions of some semilinear parabolic problems. 
Ann. Inst. H. Poincare Anal. Non Lineaire, 16:49-105, 1999. 

[16] J. P. Incardona, J. H. Lee, C. P. Robertson, K. Enga, R. P. Kapur, and H. Roelink. 
Receptor-mediated endocytosis of soluble and membrane-tethered sonic hedgehog by 
patched-1. Proc. Natl. Acad. Sci. USA, 97:12044-12049, 2000. 

[17] S. Kamin and L. A. Peletier. Singular solutions of the heat equation with absorption. 
Proc. Amer. Math. Soc, 95:205-210, 1985. 

[18] B. Kawohl. When are solutions to nonlinear elliptic boundary value problems convex? 
Comm. Partial Differential Equations, 10:1213-1225, 1985. 

[19] D. Kinderlehrer and G. Stampacchia. An Introduction to Variational Inequalities and 
Their Applications. Academic Press, New York, 1980. 

[20] A. D. Lander, W. C. Lo, Q. Nie, and F. Y. Wan. The measure of success: constraints, 
objectives, and tradeoffs in morphogen-mediated patterning. Cold Spring Harbor Per- 
spectives in Biology, l:a002022, 2009. 

[21] M. Lucia, C. B. Muratov, and M. Novaga. Linear vs. nonlinear selection for the 
propagation speed of the solutions of scalar reaction-diffusion equations invading an 
unstable equilibrium. Commun. Pure Appl. Math., 57:616-636, 2004. 

[22] A. Martinez-Arias and A. Stewart. Molecular principles of animal development. Ox- 
ford University Press, New York, 2002. 

[23] C. B. Muratov, P. V. Gordon, and S. Y. Shvartsman. Self-similar dynamics of mor- 
phogen gradients. Phys. Rev. E, 84:041916 pp. 1-4, 2011. 

[24] L. Oswald. Isolated positive singularities for a nonlinear heat equation. Houston J. 
Math., 14:543-572, 1988. 

[25] H. G. Othmer, K. Painter, D. Umulis, and C. Xue. The intersection of theory and 
application in elucidating pattern formation in developmental biology. Math. Model. 
Nat. Phenom., 4:3-82, 2009. 



16 



[26] M. H. Protter and H. F. Weinberger. Maximum principles in differential equations. 
Springer- Verlag, New York, 1984. 

[27] G. T. Reeves, C. B. Muratov, T. Schiipbach, and S. Y. Shvartsman. Quantitative 
models of developmental pattern formation. Devel. Cell, 11:289-300, 2006. 

[28] Giovanni Sansone. Equazioni Differenziali nel Campo Reale, Vol. 2. Nicola Zanichelli, 
Bologna, 1949. 2d ed. 

[29] L. Veron. A note on maximal solutions of nonlinear parabolic equations with absorp- 
tion. arXiv:0906.0669v2 [math.AP], 2011. 

[30] O. Wartlick, A. Kicheva, and M. Gonzalez-Gaitan. Morphogen gradient formation. 
Cold Spring Harbor Perspectives in Biology, l(3):a001255, 2009. 

[31] C. E. Wayne. Invariant manifolds for parabolic partial differential equations on un- 
bounded domains. Arch. Rational Mech. Anal, 138:279-306, 1997. 



17 



